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Abstract. For a two-dimensional steady supersonic Euler flow past a convex cornered wall with 
right angle, a characteristic discontinuity (vortex sheet and/or entropy wave) is generated, which 
separates the supersonic flow from the gas at rest (hence subsonic). We proved that such a transonic 
characteristic discontinuity is structurally stable under small perturbations of the upstream super- 
sonic flow in BV. The existence of a weak entropy solution and Lipschitz continuous free boundary 
(i.e. characteristic discontinuity) is established. To achieve this, the problem is formulated as a free 
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, boundary problem for a nonstrictly hyperbolic system of conservation laws; and the free boundary 

problem is then solved by analyzing nonlinear wave interactions and employing the front tracking 
method. 
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1. Introduction and Main Theorem 



We are concerned with the structural stability of transonic characteristic discontinuities in two- 
dimensional steady full compressible Euler flows, which separate supersonic flows from the static 
gases (that is, flows with zero- velocity, hence subsonic, cf. Figured]) under small perturbations in 
the space of functions of bounded variation of the upstream supersonic flow. The flow is governed 
by the two-dimensional full Euler system, consisting of the conservation laws of mass, momentum, 



00 

and energy: 



d x (pu) + d y (pv) = 0, 

4- d.,(ouv) = f). 

(1.1) 



d x (pu 2 + p) + dy(puv) = 0, 



d x (puv) + d y (pv 2 + p) = 0, 

As usual, the unknowns u = (u,v),p, and p are respectively the velocity, the pressure, and the 
density of the flow, and 

E= ~(u 2 + v 2 ) + e{p, p) 

is the total energy per unit mass with the internal energy e(p,p). Let S be the entropy. For 
polytropic gas, the constitutive relations are 

7 / S \ (7 - 1)P 

p = Kp'exp{ — J, e= 
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for some positive constants k, c u , and 7 > 1. The sonic speed is given by 




The flow is said to be supersonic (resp. subsonic) at a state point if u 2 + v 2 > c 2 (resp. u 2 + v 2 < c 2 ) 
there. It is well-known that the Euler system (jl.ip is hyperbolic for supersonic flow, and particularly 
hyperbolic in the positive x-direction if u > c; while it is of hyperbolic-elliptic composite-mixed 
type if the flow is subsonic. Hereafter, we use U = (u,v,p,p) to represent the state of the flow 
under consideration. 

An important physical case in which a characteristic discontinuity is generated is as follows: the 
characteristic discontinuity is a straight line emerging from a corner O (that is the positive rc-axis); 
the gas flow above (i.e., in {x £ R, y > 0}) is a uniform supersonic flow with the velocity (u, 0), 
pressure p, and density p + such that u > c + for the sonic speed c + > 0; below the characteristic 
discontinuity (i.e., in {x > 0, y < 0}), the gas is at rest with zero- velocity, pressure p, and density 
p~ . The question is whether such a transonic characteristic discontinuity is structurally stable 
under small perturbations of the upstream supersonic flow in the framework of two-dimensional 
steady full Euler equations, as shown in Figure [U Notice that the characteristic discontinuity is 
either a combination of a vortex sheet and an entropy wave or one of them. 

For related cases, when the flows on both sides of the characteristic discontinuity are supersonic, 
it has been shown to be structurally stable by Chen-Zhang-Zhu [3] in the framework of weak en- 
tropy solutions, and the L 1 -stability also holds as established by Chen-Kukreja [5]; when the flow 
is in an infinite duct and on both sides of the characteristic discontinuity the flows are subsonic, 
Bae [T] proved that it is stable under small perturbations of the walls of the duct. Characteristic 
discontinuities appear ubiquitously in Mach reflection and refraction/reflection of shock upon an 
interface. For such problems, Chen [6] and Chen-Fang [7] studied the stability of subsonic charac- 
teristic discontinuities; Fang- Wang- Yuan |9j showed the local stability of supersonic characteristic 
discontinuity in the framework of classical solutions. Also see Zhang [15] for supersonic potential 
flows past a convex cornered bending wall and related geometry. As far as we know, there have been 
no results available so far concerning transonic characteristic discontinuities when the supersonic 
flows are not C 1 but only belong to the space of functions of bounded variation. 

We remark that considerable progress has been made on the existence and stability of multidi- 
mensional transonic shocks in steady full Euler flows (see, for example, [H El fT2| [T3"l I14j ; also cf. 
[8]). In these papers, the smooth supersonic flow is given, and the key point is to solve a one-phase 
elliptic free boundary problem. However, in order to solve the perturbed characteristic discontinu- 
ity in this paper, the key point is to solve a hyperbolic free boundary problem in the framework of 
weak entropy solutions. 

In the following, we first formulate the aforementioned stability problem for the characteristic 
discontinuity as a free boundary problem for the Euler equations. Then, in Sections 2-5, we 
establish the existence and stability of the free boundary, by a front tracking method (cf. [21 [HI [11] ) . 

To this end, we now introduce characteristic discontinuities, a kind of discontinuities that sepa- 
rate piecewise classical/weak solutions of (jl.ip . Suppose that T is a Lipschitz curve with normal 
n = (m, n-i) in the plane, and the flows U = (u, v,p, p) on both sides of V satisfy the Euler equations 
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Supersonic flow 



Supersonic flow U = (u,v,p,p) 



Characteristic discontinuity 




Still gas U_~ = (0,0, p, p~) 



FIGURE 1. A characteristic discontinuity emerged from the corner O that separates the 
static gas with zero-velocity below from the supersonic flow above. 

(jl.ip in the classical/ weak sense. Then U is a weak solution to (jl.ip provided it satisfies (jl.ip on 
either side of T in the classical/weak sense, and the following Rankine-Hugoniot jump conditions 
hold along V: 

( , 



\pv\n\ + [pv]n 2 = 0, 
[pu 2 + p]ni + [puv]n 2 
[puv]ni + [pv 2 + p]n 2 



0, 
0, 



(1.2) 



[pu(E + E)]m 



\pv(E + Z)]n 2 = 0, 



where [•] denotes the jump of the quantity across T. Such a discontinuity T is called a characteristic 
discontinuity if the mass flux m = pu ■ n = (pu)ni + (pv)n 2 through T is zero. For a characteristic 
discontinuity, the first and fourth condition ([pu ■ n(E + £)] = 0) in (II. 2j) hold trivially, while the 
second ([upu ■ n] + [p]ni = 0) and the third ([vpu ■ n] + [p]n 2 = 0) imply [p] = 0. Thus, we see that, 
for a characteristic discontinuity, the only jump conditions should be 



[p] =0 



and 



u • n 



0. 



11-3) 



This implies that there might be jumps of the tangential velocity and the entropy (i.e., the density). 
Therefore, in general, a characteristic discontinuity in full Euler flow is either a vortex sheet or an 
entropy wave. We also note that (II. 3h implies (|1.2j) . 

Consider the Cauchy problem of the hyperbolic-elliptic composite-mixed system (II. ip : 




in x > 0, i/SR, 
U , x = 0, y > 0, 
ur, x = 0, y < 0. 



(1.4) 



The discontinuous function: 



U 



ur = (o,o,p,/r 



x > 0, y > 0, 
x > 0, y < 0, 
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with u > c + = y7p/p + is a characteristic discontinuity of when Uo = U_ + , and f/ _ is the 

state of the static gas below {x > 0, y = 0}. 

A weak entropy solution to problem (jl.4p can be defined in the standard way (cf. Definition ll.il 
below): In particular, it is defined as in (jl.6p - (|1.10p . but the domain of integration Q is replaced 
by {x > 0, y E M}, £ is replaced by {x = 0, y E E}, and the right-hand sides of (|1.7p ^ (|1.8p are 
replaced by zero. 

We note that the state of the static gas C/~ should be unchanged under the perturbation of the 
supersonic flow. This is a merit of such a transonic characteristic discontinuity, which enables us 
to reduce the above problem to an initial-free boundary problem of the hyperbolic Euler equations. 

Suppose that the characteristic discontinuity T is given by the equation: 

y = g(x) for x > 0, 

with g(0) = 0. Then 

n _ (s'(z)-i) 



The domain bounded by T and T, = {(x,y) : x = 0,y > 0} is written as O. We formulate the 
following free boundary problem of (jl.ip in Q: 

U = Uo on S, 

p = p on T, (1.5) 

v = g'{x)u on T, 

where the first is the initial data and the last two conditions on T come from (11.311. 



Definition 1.1. A pair (g,U) with y = g{x) £ Lip([0, oo);R) and U = (u,v,p,p) E L°°(fi;lR 4 ) is 
called a wea/; entropy solution to problem (|1.5p provided the following hold: 

^> U is a weak solution to (jl.ip in Q and satisfies the initial-boundary conditions in the trace 
sense: For any <f> E C~(1R 2 ), 

/ i^pudxcj) + pvd y 4>) dxdy + / pu(f)dy = 0, (1.6) 
Jn Jt, 

((pu 2 + p)d x (f) + puvd y (f)) dxdy + j (pu 2 + p)4> dy = p j (fini ds, (1-7) 
((puv)d x (fr + (pv 2 +p)d y (j)) dxdy + / (puv)4>dy = p / </>n 2 ds, (1-8) 



/ (pu(E + ^)a iE + H^ + -)^)dxdy+ [ pu(E + ^)4>dy = 0- (1.9) 
P P Jt P 

f/ satisfies the entropy inequality, i.e., the steady Clausius inequality: 

d x {puS) + dy(pvS) > 
in the sense of distribution in $7: For any (j) E Cg°(M 2 ) with <f> > 0: 

/ (puSd x 4> + pvSd y 4>) dxdy + / puS^dy < 0. (1.10) 

JS7 JE 
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We remark that, if (g, U) is a weak entropy solution to problem (jl.5p . then 

_ | U in {y > g(x), x > 0}, 

\u_ m {y < g{ x ), x > 0} 

is a weak entropy solution to problem (jl.4p . This can be checked by integration by parts in 
{x > 0, y < g(x)}; thus, we omit the details. From now on, we focus on the solution of problem 
(|1.5j) . The main result of this paper is the following. 

Theorem 1.1. There exists positive constants e and C depending only on Uj^ so that, if 

\\U -U + \\ By{r) <e, 

then problem (|1.5p has a weak entropy solution (g,U). Moreover, the solution satisfies 

(i) g G Lip([0, oo); R) with g(0) = and ||</||£oo[q )00 ) < Ce; 

(ii) There exists U_ G M 4 so that 

U-U G C([0, oo); L 1 ^), oo)), ||(C7 - OLvflgM.oo)) < Cs - 

Remark 1.1. We note that ||[/o — IZ + || BV (£) — £ im P nes t na t lirm y _ ! . OC) ([/o — LL + ){y) exists. Then 
there exists U_ G M 4 as claimed in Theorem 11.11 so that 

lim U (y) = C/ , 

and 

\Uo -H + \ < e- 

To prove Theorem ll.H we establish the compactness and convergence of approximate free bound- 
aries to the free boundary of the exact solution in supersonic-subsonic flows in the framework of 
front tracking method, while some other essential tools/notions of the front tracking method are 
extended, modified, and further clarified working in the presence of the free boundary such as a 
generation of fronts to control the finiteness of physical fronts and the errors from approximate 
Riemann solvers for the nonstrictly hyperbolic free boundary problem. For this, two new nonlinear 
Riemann problems are involved: One is the Riemann problem at the convex corner connected with 
the still gas state (subsonic state); and the other is the Riemann problem determining the evolution 
of the free boundary, for which we establish the boundedness of the key reflection coefficient of the 
reflected wave into the supersonic region after the interaction of the incident wave with the free 
boundary. To achieve the compactness, we have to identify the right scales and global weights 
to control the Glimm functional to make it monotonically decrease in the flow direction, while 
preserving the overall structural stability of the characteristic boundary as the hyperbolic region 
evolves in complicated ways under any small BV perturbation yet the subsonic state remains stable 
beneath the free boundary. 

We also remark in passing that, as an example of one-phase hyperbolic free boundary problems 
for nonstrictly hyperbolic systems, we deal with the problem in the physical space, the Euler 
coordinates throughout this paper. This represents a first example of an approach to apply the 
front-tracking method to study structural stability of interfaces between different mediums, one 
of them is subsonic. Our approach offers further opportunities to initiate the study of vortex 
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sheets/entropy waves in the space of bounded variation in nozzles, jets, etc. for mixed-type flows, 
transonic flows. In a forthcoming paper, we will deal with this problem and related L 1 -stability in 
a different approach. 

The rest of this paper is devoted to establishing Theorem 11.11 We will mainly employ a version 
of the front tracking method introduced in Holden-Risebro [11] for convenience to deal with the 
problem. Thus, in Section 2, we review some facts concerning the solvability of various Riemann 
problems for the steady Euler equations, and present some essential interaction estimates. It 
manifests clearly in the simplest case how such a hyperbolic free boundary problem can be solved. 
Then, in Section 3, we construct approximate solutions by the front tracking algorithm. The key 
point is to show such an approximate solution can be established for x £ [0, oo) by constructing 
a Glimm functional. Then, in Section 4, with the uniform BV estimate of approximate solutions 
obtained from the Glimm functional, we establish the compactness of the family of approximate 
solutions and that the limit is actually an entropy solution. Finally, we discuss the asymptotic 
behavior of the weak entropy solutions as x — > oo in Section 5. 

2. Riemann problems and interaction estimates 

In this section we first review certain basic properties of the steady hyperbolic Euler equations 
(jl.ip that are used later for self-containedness (cf. Chen-Zhang-Zhu [3, pp. 1665-1670]). Then we 
show the solvability of "free boundary" Riemann problem and interaction estimate between weak 
waves and the free boundary, which are the new ingredients in this paper. 

2.1. Euler Equations. As in [3], we write the Euler equations (|1.1|) in the form 

d x W(U) + d y H(U) = 0, U=(u,v,p,p), (2.1) 

where 

9 9 

W(U) = {pu,pu 2 +p,puv,pu(- — H — )) T 

(7 - l)p 2 



and 



H(U) = (pv, puv, pv 2 +p,pv( + — - — )) T . 



The eigenvalues A of this system are determined by det{W jjW{U) — VjjH(U)) = 0, or explicitly, 

(v - Xu) 2 ((v - Xuf - c 2 (l + A 2 )) = 0. 
Thus, if u > c, we have four real eigenvalues: 

u z — c z u 
The associated linearly independent right-eigenvectors are 

r j = K j (-X j ,l,p(X j u-v),^^^-) T , i = 1,4; (2.3) 

r 2 = (n, U ,0,0) T , r 3 = (0,0,0, /3 ) T , (2.4) 

where Kj are renormalized factors so that r,- • WjjXj(U) = 1 since the j-th characteristic fields 
are genuinely nonlinear, j = 1,4. While the second and third characteristic fields are linearly 
degenerate: • VfjXj(U) = 0, j = 2, 3. Although the steady Euler system is not strictly hyperbolic, 
we can still employ the general ideas presented in [8j [TT] to treat related Riemann and Cauchy 
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problems. The only difference is that, although the characteristic discontinuity has only one front 
in physical space (since two of the four characteristic eigenvalues coincide) , we need two independent 
parameters (one corresponds to A2 for the vortex sheet, and the other to A3 for the entropy wave) 
to represent its strength. 

At the unperturbed reference state U_ + = (u, 0,p, p + ), we easily see that 

AiO/ + ) < A 2 (£/+) = = X 3 (U + ) < X 4 (U + ) = -AiO/ + ). 

Also, Lemma 2.3 in [3] indicates that the re-normalization factors Kj(U),j = 1,4, are positive in a 
small neighborhood of U_ + . 

2.2. Wave Curves in the Phase Space. As shown in [3j, at each state Uq = (uq, vq,pq, po) with 
uq > Co in the phase space, there are four curves in a neighborhood of Uq: 

Vortex sheet curve C^fTo) : U = (uoe a2 ,voe a2 ,po, pq). 

These are the states U that can be connected to Uq by a vortex sheet with slope ^ and 
strength 0-2 £ R; 
Entropy wave curve C 3 {Uq) : U = (uq, vq,pq, poe as ). 

These are the states U that can be connected to Uq by an entropy wave with slope — 
and strength 03 € R. 
<C> Rarefaction wave curve Rj(Uo): 

dp = c 2 dp, du = —Xjdv, p(XjU — v)dv = dp for p < pq, u > c, j = 1,4. 

These are the states U that can be connected to Uq from the lower by a rarefaction wave 
of the j-th family; 
Shock wave curve Sj(Uo): 
c 2 

\p\ = r[/»]i M = -8j[v],po(sjU - v )[v] = [p] for p > p ,u > c, j = 1,4. 

These are the states U that can be connected to Uq from the lower by a shock wave of 
the j-th family, with the slope of the discontinuity to be 

_ ^o^o + (-Ipcy^p + v'q ~ c 2 . _ 

s i ~ ,,2 ^2 ' — ' ' 

u — c 

where c = ^4 and b = ^t- — . 

PO" 2 2 po 

One can also parameterize Rj(Uo) and Sj(Uo) (j = 1,4) so that there is a curve given by a C 2 
map aj h-> <&j(otj; Uq) in a neighborhood of Uq, with ay > being the part of Rj(Uq), and ay < 
the part of Sj(Uo), and 

$,(0; Uq) = Uq, d a .^(0; Uq) = rjiUo). (2.5) 

We can also write the curve Cj(Uo) (j = 2,3) as aj 1— > $j(aj;Uo) which is still C 2 so that (|2.5p 
hold for j = 2, 3. Since {r 3 -(Z7o)}| =1 are linearly independent, such curves consist locally a (curved) 
coordinate system in a neighborhood of Uq. This guarantees the solvability of the Riemann problems 
stated below. 

For simplicity, we set 

$(a 4 > 0-3,0*2 ,0!i;[/q) = $4(04; $3(0-3; $2(0-2; $i(oi;C/ )))). (2.6) 
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Then 

$(0,0,0,0;£7 ) = *7o, d a ^(0,0,0,0;U )=r j (U ), j = 1,2,3,4. (2.7) 

2.3. Standard Riemann Problem. We now consider the standard Riemann problem, that is, 
system (jl.ip with the piecewise constant (supersonic) initial data 

{U + , y > yo, 
(2.8) 
U , y < 2/o, 

where U + and U~ are the constant states which are regarded as the above state and below state 
with respect to the line y = yo, respectively. 

Lemma 2.1 (Lemma 2.2 in [3]). There exists e > such that, for any states U~ and U + lie in 
the ball € {Uq) C M 4 with radius e and center Uq, the above Riemann problem admits a unique ad- 
missible solution consisting of four elementary waves. In addition, the state U + can be represented 
by 

£7+ = <&(a4,a 3 ,a 2 ,ai;U~). (2.9) 

It is noted (cf. Lemma 4.1 in [3]) that one can use the parameters ctj,j = 1, ... ,4, to bound 
\U + — U~\: There is a constant B depending continuously on Uq and e so that, for J7 ± connected 
by (USD, 

1 4 4 

-Y,\^\<\u + -u~\<b^2\^\- 

3=1 3=1 

For later applications, it is also important to express the Riemann solver from the upper state 
U + to the lower state U~, rather than the usual way given above. For U + = &j(otj;U~), we 
may have a C 2 -map U~ = * j (a j ;U + ) with Vj(0;U) = U and d aj ^j(0;U) = -Tj(U). Thus, for 
U + = <I > (q4, 03, Q2, ai; U~), we may express U~ in terms of U + by 

LT = *(ai, q 2 , a 3 , a 4 ; t/ + ) = *i(ai; ^ 2 (a25 ^3(03; ^4(04; ^ + ))))- 
Then ^(0,0,0,0; 17) = U and 5 Oj .*(0, 0,0,0; 17) = -Tj(U). 

2.4. Free Boundary Riemann Problem. We now consider the following Riemann problem of 
(jl.ip involving a free boundary — a characteristic discontinuity. The initial data is a constant state 
U = U + given on the positive y-axis, and the free boundary is a straight line y = kx with k G R 
to be solved. The boundary conditions on the free boundary are p = p and fc = -. Since the free 
boundary - characteristic discontinuity — is of the second/third characteristic family, the Riemann 
solver should contain only one 4- wave with parameter 0:4 and a middle constant state U*; see Figure 
2 below. 

Lemma 2.2. There exists e > so that, for U + G O e (U_ + ), there is only one admissible solution 
consisting of a J^-wave that solves the above free boundary Riemann problem. The middle state U* 
can be represented by U* = ^4(04; U + ), and the free boundary is determined by k = There also 
holds 

a 4 = K 1 (p + -p) + M 1 \U + -U + \ 2 , \k\ < K[\U + - U + \, (2.10) 
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Characteristic discontinuity 



FIGURE 2 . A Riemann problem with a free boundary that is a characteristic discontinuity. 
with the constants K\ , K[ > and a bounded quantity Mi only depending continuously on U_ + and 



Proof. 1. We write to denote the k-th argument of the vector U, k = 1, . . . , 4. Consider the 
function: 

L(a,U + ) = (y 4 (a;U + ))^ -p = (V 4 (a;U + ) - * 4 (0;I/ + )) (3) , 

for which L(0;U + ) = 0. Then 

d a L(0;U + ) = -(r 4 (C/ + )) (3) = -(k aP uX a )\ e + < 0. 

From the implicit function theorem, we infer that a can be viewed as a function of U + G O t (U_ + ) 
for suitably small e > 0. In particular, a(U_ + ) = 0. This completes the existence proof. 

2. Since Vt/tf 4 (0; tf) = h, duL{0;U + ) = (0,0, 1,0). Then 

(0,0,1,0) 



Vua{ir 

Thus, by the Taylor expansion, we conclude 



(k 4 ^A 4 )| 



u+ 



where Ki 



i 



ll (k4PmA4)| [7 

3. From the above, we have 
Then we have 



a = K!(p + -p) + M l \U + -U + \ 2 , 
> 0, and Mi is a constant depending continuously and only on U_ + and e. 



\U* - U + \ < B\a\ < B'\U + - U + \. 



\U* - U+\ < B"\U + - U + \ 

for some constant B" > 0. Hence, regarding ^ as a function of U and by the mean value theorem, 
we have 

| — | < C\U* - U + \ < K[\U + -U + \ 
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as desired. □ 

2.5. Approximate Riemann Solver. The front tracking method involves approximating the 
rarefaction waves appeared in the Riemann problems or (free) boundary Riemann problems by 
several artificial discontinuities separating piecewise constant states. 

Suppose that U + = $(0:4, 03, Q2, a\\ U~) gives the solution to the standard Riemann problem 
(j!T8]) . with middle states U 1 = $i(ai;U~) and U 2 = ^ 4l (a 4l ;U + ). For any 5 > 0, we define a 
5-approximate solution II s to the Riemann problem as follows: 

• If a\ > 0, then the 1-wave is a rarefaction wave that requires modification as follows. Set 
v be the closest integer to ^ (that is, v G Z and ^ — \ < v < + 5), as well as 
#1,0 = U~, U ljU = U 1 , and Ui yk = $i(^ai; for k G {1, . . . , ^ - 1}. Then, in the 

wedge {(x,y) : x > 0, y < A*x}, we define 



U 5 



U-, y<Ai(C7-K 

^i,fc, Ai(£7i jfc _i)x < y < Xi(U 1>k )x, k = 1, . . . , v - 1, (2.11) 
U 1 , Xi(Ui tV -i)x < y < A*x. 

Here A* is a constant chosen so that sup Ue0 ^ u +^ X\ < A* < inf^Q m+) A2, which exists 
when e is small. 

Then the rarefaction wave is replaced by "step" functions with width (strength) Si, and 
the discontinuity between two steps moves with the characteristic speed of the lower state. 
If a 1 < 0, then the 1-wave is a shock, and no change is necessary. In the wedge {(x,y) : 
x > 0, y < A*x}, we define 



U 6 



U , y < six, 
U 1 , s\x < y < A*x, 

where s± is the speed of the shock front. 

• For a2,«3, there is always no change. 

• Similar to the case of the 1-wave, we can define U s in {x > 0,y > — A*x} by consid- 
ering whether the 4-waves is a rarefaction wave (with modification) or a shock (without 
modification). 

2.6. Interaction of Weak Waves. The following weak wave interaction estimate is classical; see 
Lemma 3.2 in [H p. 1670]. 

Lemma 2.3. Suppose that U + ,U m , and U~ are three states in a small neighborhood of Uq with 
U+ = ^,a 3 ,a 2 , ai ;U m ), U m = $(/3 4 , fo, U~), and U+ = $(74, 73, 72, 7i5 U~). Then 

l j = a J +(3 j + 0(l)A(a,(3), (2.12) 
where A(a,/3) = |a 4 |(|/3i| + |/3 2 | + + (|a 2 | + \a s \)\Pi\ + Ej=i,4 A j(«> Z 3 )* with 



A,-(a,/3) 



0, aj > 0, p* > 



I ay 1 1 /3j |, otherwise. 
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2.7. Interaction of Weak Wave and Free Boundary. We now consider the change of strength 
when a weak wave interacts with the free boundary (see Figure 3). It is only possible that a 
weak 1-wave at impinges on the characteristic discontinuity S , and resulting a reflected 4- wave 
with parameter 04, and the characteristic discontinuity itself is also deflected to a new direction, 
denoted to be S r . We note that both U r and S r can be solved by the free boundary Riemann 
problem with initial data U m . 




FIGURE 3. A 1-wave ol\ is reflected by the characteristic discontinuity S , resulting in a 
reflected 4- wave «4 and deflected characteristic discontinuity S r . 

Lemma 2.4. Suppose that U , U m , and U r are three states in O e (U_ + ) for sufficiently small e, with 
U m = $i(ai; U l ) = $4(04; U r ). Then 

a 4 = -K 2 a 1 + M 2 \ai\ 2 , (2.13) 

with the constant K 2 > and the quantity M 2 bounded in O e (U_ + ). Furthermore, for U = 
(ui,vi,pi, pi), \K 2 \ > 1, \K 2 \ < 1, and \K 2 \ = 1 when v\ <0, vi > 0, and v\ = 0, respectively. 

Proof. 1. We have U m = $i(a; C/') and U r = ^4(f3; U m ). Consider the following function: 

L{fi,a) := (* 4 (A^l(«;^))-^) (3) - 

Then L(0, 0) = 0, and dpL(0,0) = —(r^U 1 ))^ < 0. By the implicit function theorem, there exists 
a function j3 = /3(a) so that L(f3(a),a) = for small a. We see /3(0) = 0. 

2. We calculate d a L(0,0) = (r^C/'))^ < 0. Thus, ^ = -if 2 := ?^t|%S < 0. Therefore the 
equality in (I2.13P follows from Taylor expansion. 

3. The coefficient 

K (n(^ 3 ) _ ^-Ai(^) 

2 - (r 4 (C/0)( 3 ) ^ + A 4 (C/0 
and, for any state U = (u,v,p,p) G O e (^ + ), there holds Xi(U) < X 2 ^(U) = ^ < ^(U). Using 
these two facts with the expressions for X\(U) and \4(U) given in (|2.2p . it follows that \K 2 \ < 
1, > 1, and \K 2 \ = 1 when v\ > 0, V\ < 0, and v\ = 0, respectively. □ 



12 



GUI-QIANG CHEN, VAIBHAV KUKREJA, AND HAIRONG YUAN 



3. Construction of Approximate Solutions and Uniform Estimates 



In this section we adopt the front tracking method in Holden-Risebro [11 J to construct a family 
of approximate solutions {(g s , U s )}s>o of the problem (|1.5|) and present some uniform estimates 
independent of 5, which is necessary for a compactness argument in §4 to show the existence of a 
weak entropy solution to (jl.5p . 

3.1. Construction of Approximate Solutions. For any given 5 > 0, we now describe the 
construction of an approximate solution (g s , II s ) to the free boundary problem (jl.5p . 

We first approximate the initial data Uo(y) by a piecewise constant function Uq{u) as done in 
the study of the Cauchy problem. We require that 



We solve the Riemann problems with initial data on {x = 0, y > 0} and a free boundary Riemann 
problem at the corner (0,0), and then approximate rarefaction waves as carried out in §2.5 with 
parameter 5 to obtain new discontinuities. Note the resulting (approximate) solution is piecewise 
constant. 

Then we need do nothing until as x increases to some value x = t, where 

(i) either two fronts interact; 

(ii) or there is a weak 1-wave that interacts the free boundary (it is obtained by solving the 
free boundary Riemann problem before) from above. 

As noted in [2], by adjusting the slopes of the discontinuities, we can assume that, at each 
{x = t}, only one of the two cases above happens. This is harmless since the error can be made to 
be arbitrarily small. 

For case (i), as mentioned above, by adjusting the slopes of these discontinuities (with arbitrarily 
small error), we may assume that only two discontinuities collide. Suppose that the lower disconti- 
nuity is of r-family and has a parameter a with the lower (constant) state U l and upper (constant) 
state U m , the upper discontinuity is of s-family and has a parameter j3 with the lower (constant) 
state U rn and upper (constant) state U r ', and they collide at the point (r, rj). Then, as before, we 
solve a Riemann problem at (r, rj) with the lower state U l and upper state U r , by applying the 
approximate Riemann solver to obtain new discontinuities. 

For case (ii), we may still assume only one discontinuity collides with the free boundary (tran- 
sonic characteristic discontinuity) . Then we solve a wave reflection-deflection problem with a 1-wave 
reflected by the free boundary, obtaining a reflected 4-wave and a deflected characteristic disconti- 
nuity (see Figure 2). If the reflected 4-wave is a rarefaction wave, by approximating the rarefaction 
wave, we obtain again the approximate solver containing new discontinuities. 

Continuing this procedure and, in some cases, removing certain quite weak fronts (cf. £13.4.21 
below for details), we obtain an approximate solution (g s , II s ). 

Remark 3.1. To ensure that the above procedure works to construct an approximate solution for 
all x £ [0, oo), we need to show that, for any < x < oo, 
• The total variation is small: T.V.(U s (x, •)) < Ce; 





By Remark 11.11 we may also assume that, for each 5 > 0, there holds U[ 



(y) = Ho for lar g e V- 
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• An L°°-bound: The solution still lies in a small neighborhood of U_ + ; 

• Given any finite T > 0, there happens only a finite number of collisions/reflections for 
{0 < x < T}. 

The first two are necessary so that we can actually solve the standard or free boundary Riemann 
problem. Here C is a universal constant independent of e and x > 0. The third one guarantees 
that the global approximate solutions defined up to any x > can be actually obtained. 
In the following three subsections, we deal with these three issues. 

3.2. Bounds of Total Variation. We now establish the bounds of total variation of the approx- 
imate solutions U s (x,y). 

3.2.1. Glimm Functional. We introduce the following version of Glimm functional 

G{x) = V{x) + kQ(x), (3.2) 

where k > is a large constant to be chosen. The terms V and Q are explained below. By the 
properties of the approximate Riemann solver, T.V.(U s (x, •)) is equivalent to V(x). Then it suffices 
to prove 

V{x) < C e (3.3) 

for a constant Co depending only on U_ + . Recall here e = \\Uq — fZ + || BV ([o oo)) measures the strength 
of the perturbation of initial data. 

For a weak wave/discontinuity a of i a -family, we define its weighted strength as 

k+a if a G T+ and i n = 1, 

(3.4) 

a if a £ Tt and i a = 2, 3, 4, 

where k + > \K%\ for the coefficient K<i appeared in Lemma 12.41 and we use Tt to denote the set of 
weak waves/discontinuities (not including the free boundary) that cross the line {x = t}. 

• The weighted strength term V(t). We define the total (weighted) strengths of weak waves/discontinuities 
at x = t as 

V(t)=J2\ba\- (3.5) 

• The interaction potential term Q(t). The interaction potential term we use here is the same 
one as introduced by Glimm [10], that is: 

Q{t) = Yl \ bab ^ ( 3 - 6 ) 

where A(t) is the approaching set defined by pairs (b a ,bp) so that, for x = t, the waves/discontinuity 
with strength b a lies in the lower side of the waves/discontinuity with strength bp, and b a is of 
family i a and bp is of family ig, where i a > ip, or both are of the same family but at least one of 
them is a shock. Note we do not consider the free boundary as a wave/discontinuity in this paper. 

As shown by Lemma 6.2 in |llj . at x = r, if two discontinuities of strengths b a and bp collide, 
then we have 

Q(r+) - Q(r-) = ~\b a bp\, (3.7) 
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provided that 

V{t-)<h:=\o{1). (3.8) 
It is here one needs Lemma 12.31 If no discontinuities collide at x = t, then Q(r+) = Q(t—). 

3.2.2. Non-increasing of the Glimm Functional. We now show the bounds of total variation by 
proving that the Glimm functional G(x) is non-increasing for x. There are the following three 
cases. 

(i) Collision of discontinuities. For x = r where two discontinuities b a and bp collide, there is 
no other wave interaction and reflection upon the free boundary as we assumed. Therefore, 
the decreasing of G{r) is classical. By Lemma |2.3^ we have 

G(t+)-G(t-) = ( V (t+)-V(t-))+k(Q(t+)-Q(t-)) 

< M\b a bp\ + K(--\b a b p \)<0, 

if we choose k > 2M sufficiently large. Note that 0(1) does not depend on the approxima- 
tion parameter 8. 

(ii) Weak 1-wave interacts with the free boundary. For weak wave a\ of 1-family 
interacts with the free boundary from above, resulting in a reflected 4-wave a±. By Lemma 
12.41 we have 

G(t+)-G(t-) = ( V (t+)-V(t-))+k(Q(t+)-Q(t-)) 

< \b ai \ - \b ai \ + Kfi\b a4 \ 

< {{nn + l){-K 2 + M 2l i)-k + )\a l \ <0 

if we choose k + sufficiently large (independent of 5). 

(iii) Other situation. If, for x = r, no collision or reflection upon the free boundary happens, 
then we still have G{t+) = G{t—). 

In the above, we have determined k and fc+ independent of 5, and proved that, for any x = r > 0, 
there holds G(r+) < G(r-), provided (JHSJ) holds. 

3.2.3. Boundedness of Total Variation. The bound V(r) < Cq£ then follows from an induction 
argument as shown in p. 217] for the proof of Lemma 6.3 there, provided that e is small. 

We first set < t% < ... < < ... as the sequence so that, for x = r^, either collision or 
reflection upon the free boundary occurs, and set Vk, Gk the value of V(tj.— ), Gfa— ) respectively. 

We know that there exists a constant C% independent of 5 > so that V(r) < CiT.V.(U <5 (r, •)) 
for all x > 0. Note here the choice of weight fc+ is in essence only determined by U + . Define 

C = C 1 + kCI 

We choose positive e < 1 small so that 

Cie + K(Cief < n, C 2 C e < e. 

Here e is the value so that the Riemann problems or the free boundary Riemann problems can be 
solved when the Riemann data are in O e (U_ + ), and C 2 is the constant depending only on U_ + so 
that T.V.{U 5 (x, •)) < C 2 V{x) for any x > 0. 
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By assumption on the initial data, we have T.V.(C/q ) < e. Thus, by a property of the Riemann 
problem, we may have 

V\ < C\E < min{Co£,/i}, 

and furthermore, 

Gi < Vi + fcV\ < Cie + n{Cie) 2 < min{C e,M}- 
Suppose that, for n < k, we have proved 

V n < min{C e,^}. 

Then, by decreasing of the Glimm functional, we have proved that there holds 

Vfc+i < Gk+i < Gf. < . . . < G\. 

This shows 

V n < min{Co£, fi} for all n. 

If we further choose e small so that CqE < fi, we obtain the bound V(r) < Cq£ as desired. This 
again implies the uniform estimate: 

T.Y.(U s (x,-))<C 2 C e. (3.9) 

3.3. L°°-Estimate of {U s } and Lipschitz Estimate of {g s }. The fact that {U s } 5>0 is uni- 
formly bounded follows directly. For each x, the solution II s (x,y) is just the constant state U_ 
for sufficiently large y, by the finiteness of propagation speed and the fact that the initial data 
^o(y) ~~ ^ f° r V ~^ 00 • Since we have proved T.V.(U s (t, •)) < C2C0S for any t > 0, then, by 
definition of the total variation, we conclude 

\\U s {x,-) -U + \\ LX <C 2 C e (3.10) 

for some new constant C 2 - 

Estimate (|3. 10[) implies the following uniform estimate on the free boundary that is given by the 
equation y = g s (x): 

\\(g 5 y\\ Loc <C 3 e. (3.11) 

with a constant C3 depending only on U_ + . In particular, by construction, for fixed 5 > 0, g s is 
a piecewise linear (affine) function, and except for countable points {r^}, it is differentiable, with 
(g s )'(x) = V s[ x ' 9 s\ X y, ■ Thus, by the mean value theorem, 

\{g S )'(x)\ < C'\U 5 (x,g\x))-U + \ < C"||[/ 5 -C/+|| Loo < C'C 2 C e, 

where the constant C depends only on U_ + . 

3.4. Finiteness of Collisions and Reflections. To show that the numbers of fronts/discontinuities 
and collisions/reflections do not approach infinity in {0 < x < r} for any finite r > 0, the basic idea 
presented in [1TJ for the Cauchy problem works well, but we have to consider additional issues such 
as the reflections off the free boundary and the fact that the Euler system is not strictly hyperbolic 
in the argument. For completeness, we give the proof below, which closely follows that in [IT] . 
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3.4.1. Generation of Fronts and Modified Construction of Approximate Solutions. Firstly, we define 
the notion of generation of a front. We set that each initial front starting at x = belongs to the 
first generation. Take two first-generation fronts of families d and h, respectively, that collide. The 
resulting fronts of families d and h belong to the first generation, while all the remaining fronts 
resulting from the collision are called second-generation fronts. Generally, if a front of family d and 
generation m interacts with a front of family h and generation n, the resulting front of families d 
and h are still of generation m and n, respectively, while the remaining fronts resulting from this 
collision are given generation n + m. The fronts of 4-family resulting from reflection of a front a of 
1-family off the free boundary has the same generation of the front a. The point of this notion is 
that the fronts of high generation are quite weak. 

Given the approximation parameter 5 > 0, we remove all fronts with generation higher than N, 
with 

N=[hnKr(S)] (3.12) 

in our construction of approximate solution (g s ,U 5 ). Here [z] denotes the integer larger than but 
closest to z and, following the notations in p. 218], we set 

r = TM=|>|<VW, K = ^±- 

with £o sufficiently small and fixed, and taking later e < £q, so that T < ™. 

More precisely, if two fronts of generation n and m collide, at most two waves will retain their 
generation. If n + m > N, then the remaining waves will be removed; however, if n + m < N, 
we use the original (approximate) solution. When we remove the fronts, we let the function II s be 
equal to the value that has to be the lower of the removed fronts, provided that the removed fronts 
are not the upmost fronts in the solution of the Riemann problem. If the upmost are removed, 
then II s is set equal to the value immediately to the upper of the removed fronts. 

We remark that this process of removing (very) weak waves in approximate Riemann solver in 
our construction of approximate solutions will not influence the uniform estimates we obtained in 
§3.2-3.3. In particular, we still have T < -X?. 



3.4.2. Finiteness of Fronts and Collisions. We will show that there exists only a finite number of 
fronts of generation less than or equal to N and that, for a fixed 5, there is only a finite number of 
collisions / reflections. 

For this, as we know that T < jg, then the strength of each individual front is bounded by ™. 
For later reference, we also note that, by f|3. 12[) . 

(AKT) N+l < 5. (3.13) 

First we consider the number of fronts of first generation. This number can increase when the 
first-generation rarefaction fronts split into several rarefaction fronts. By the term rarefaction front 
we mean a front approximating a rarefaction wave. Note that, by the construction of the approxi- 
mate Riemann problem, the strength of each split rarefaction front is at least |<5. Given that T is 
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uniformly bounded, we find 

4T 

(Number of first — generation fronts) < (Number of initial fronts) H -. (3-14) 

3o 

Thus, the number of first-generation fronts is finite. This also means that there will be only a 
finite number of collisions/reflections between first-generation fronts and free boundary. To see this, 
note first having the assumption of strict hyperbolicity would have implied that each wave family 
will have speeds that are distinct. However, we see that, although the Euler system is not strictly 
hyperbolic, the multiplicity of the eigenvalues is constant for the states U near the background 
state U + . That is, \i(U) < X 2 (U) = X 3 (U) < A 4 (C7) for any state U G O e (U + ), and hence the 
eigenvalues are separable in the same way for any state U. 

Hence, we can still conclude that each first-generation front will remain in a wedge in the (x, y)— 
plane determined by the slowest and fastest speeds of that family. Eventually, all first-generation 
fronts will have interacted at most finite times, and we can also conclude that there can be only a 
finite number of collisions between first-generation fronts and free boundary globally, since once a 
front is reflected, it will never meet the free boundary again. 

Assuming now that, for some m > 1, there will be only a finite number of fronts of generation 
i, for all i < m, and that there will only be a finite number of interactions between the fronts and 
fronts reflection off free boundary of generation less than m. Then, in analogy to (|3,14p . we find 

Number of m-th generation fronts 

4T 

< 2 x (Number of j-th and i-th-generation fronts; i + j = m) H < oo. (3.15) 

3d 

Consequently, the number of fronts of generation less than or equal to m is finite. We can now 
repeat the arguments above showing that there is only a finite number of collisions between the 
first-generation fronts (and reflections off free boundary), just replacing "first generation" by "of 
generation less than or equal to to" and show that there is only a finite number of collisions 
producing the fronts of generation of m + 1. Thus, we can conclude that there is only a finite 
number of fronts of generation less than N + 1, and that these interact (reflect off free boundary) 
only a finite number of times. 

4. Convergence and Existence of Weak Entropy Solutions 

In this section we show the strong convergence of a subsequence of the approximate solutions to 
a weak entropy solution of problem (jl.5p . 

4.1. Compactness. We first show there exists a subsequence of approximate solutions {(g s , U s )}s>o 
that converges to some (g, U) almost everywhere. In §4.2, we show that (g, U) is actually a weak 
entropy solution to problem ([1.5P . 

4.1.1. Compactness of {g 6 }- We first show the compactness of the approximate free boundary 
{fl r<5 }(5>0- More explicitly, we have 

Lemma 4.1. Let g s (x) be the free boundary for the approximate solution U s (x). Then there is a 
subsequence 5j — > so that g Sj (x) — > g(x) uniformly in any compact set. Furthermore, the limit 
g{x) is Lipschitz continuous: \g(xi) — g(x2)\ < C^e\x\ — for some constant C3. 
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Proof. By (|3.1ip . that is, IIG^yilL^Qo.oo)) — C%e and g s (0) = 0, we see that, for fixed T > 0, 
the family {g s } is uniformly bounded and equicontinuous on [0, T]. Then, by the Arzela-Ascoli 
compactness criterion, there is a subsequence Sj — > so that g s i =4 g uniformly for some g in [0, T] 
and one easily proves that \g(xi) — g(x2)\ < C^e\x\ — x 2 \ for x\,x 2 G [0,T]. By taking a diagonal 
subsequence for 2T, 3T, . . ., we can prove that g is defined for x G [0, oo) and g Sj — > g uniformly in 
any compact subset of [0, oo), and \g(xi) — g(x%)\ < C3e|xi — x 2 \ for any finite x% and x 2 . □ 

4.1.2. Compactness of {II s }. We use the following compactness lemma, which is a modification of 
Theorem A. 8 in [llj. 

Lemma 4.2. Let {u v : [0, 00) x [0, 00) — > M. }^ be a family of functions such that, for each positive 
T, 

(a) \ur){x, 0)\ < Ct for (x,8) £ [0,T] x [0, 00) u>i£/i a constant Ct independent ofr]; 

(b) For all t £ [0,T], there holds 



sup / \u ri (x,9 + ^) - u v (x,9)\ d9 < vb,t(\p\), 
\Z\< p Jb 

for a modulus of continuity v and all compact B C [0, 00) (here u v (x,t) is extended to be 
zero for x ^ [0, 00) ); 
(c) Furthermore, for any R > 0, for s and t in [0,T], there holds 

|iijj(t, 9) — u v (s, 9)\d9 < ujT{\t — s\) as 77 — > 0, 

for some modulus of continuity lot- 
Then there exists a sequence rjj — > such that, for each x G [0,T] ; the function u^^x) converges to 
a function u{x) in L 1 ([0,oo)). The convergence is in the topology o/C([0,T];L 1 [0,oo)). 

For any T > 0, note that U s (x,y) is defined for < x < T and g s (x) < y < 00. By introducing 
9 = y — g S (x), we may regard U s as a function of 9 G [0,oo) and x G [0, T] by defining 

U s (x,9) = U s {x,9 + g 5 (x)) 



^ = \\U s \\ Loa , and T.Y.(U s )(x,-) = T.Y.(U s )(x,-). Then, 



to apply Lemma I4T21 Obviously U 

■■ ■■ — 

by (|3.10p . we see immediately that (a) is valid for {U s } s>0 - 

Using the boundedness of L°° norm and total variation of U s (cf. (|3.9fl ) . the verification of (b) is 
elementary. Without loss of generality, we assume £ > 0. Then by monotone convergence theorem, 

•(fc+i)« 



/ \U s {x,9 + i) - U s (x,9)\d9 = V / \U 5 {x,9 + £) - U 5 (x,9)\d9 
Jr+ k=0 Jkt 

£ OO 

= / V + + 1)£) -U s (x,z + k£)\dz 



k=0 

<(T.V.U\x r )M\<(C 2 C eM\. 

The verification of (c) is also not difficult. For < s < t < T, we will prove that 

-R 

\U s (t, 9) - II s (s, 9)\d9< C(t - s), (4.1) 
for any R > and a constant C independent of 5,t, and s. 
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To this end, for given approximate solution U s , suppose the "collision times" are 

o < n < . . . < r fc < . . . . 

Then, for x G (rj,Tj + i), nothing happens on the (approximate) free boundary, and then we may 
ignore the free boundary and write II s (x,y) in the form 



U s (x, y) = Y,Vi+x ~ Ul)H(y - y{{x)) + U{, 



(4.2) 



k=l 



with H(-) the Heaviside step function (whose value is for the negative argument and is 1 for the 
positive argument). Here we have assumed that, for x E (rj,rj+i), there are iVj discontinuities with 
equation y = x' l k (x) (from the lower to upper as k = 1, ... , iVj), and the state in the lower side of 
{y = x l k } is U k . Prom ^3.4.2| we know that < oo. 

With the above expression, for Tj < s < t < r i+1 , we have 



\U\t,6)- U d {s,6)\d9 



de 



dg S (r) 






) 


dr 


+ 


dr 





< / / f^\Ut +1 -Ul\\H'((g 5 (r)+e)-yl(r)) 

jm+ Js 

<(L + C 3 e) / £1^+1-^1 / \H'((9 S (r) + 6)- y*(r))| dd dr 
Js r^i Jr+ 



drd8 



k=l 



(L + C 3 e) 



ft 

/ El^+i-^l dr 
Js k=i 

< (L + C7 3 e)T.V.([/' 5 (r i +,-))(t-s) < {L + C 3 e)C 2 C e(t - s). 
Here we have set 



(4.3) 



L= sup (|A 1 (C/)|,|A 2 , 3 (C/)|,|A 4 (C/)|) 
ueo e (u+) 



to be the maximal characteristic speed, and used the fact that 

dr 



used to control I ty- 1 



dx 



(4.4) 

< L. Estimate (j3.11[) is also 



We note f|4.3f) also holds for s = t% and/or t = r^+i. Then, for s G (r^r^+i) and i G ( r j) r j+i) 
with i < j, using ()4.3|) repeatedly in the intervals (s, 73+1), (r^+i, ri+2), • • • , (r^-i, rj), and (rj,t), we 
obtain flU) with C = (L + C 3 £)C2C £. 

Therefore, by Lemma 14.21 we can find a subsequence {f/^} that converges to some U under 
the metric of C([0, T]; L 1 ([0, 00))). In addition, upon at most a further subsequence, g Sj — > g. 
Now set U(x,y) = U(x,y — g(x)), which is defined in the domain Q = {x > 0,y > g(x)}, with 
D = {y = g(x)} being the lateral (free) boundary. In §4.2, we show that (g, U) is actually a weak 
entropy solution of problem (jl.5p . In the following, for simplification, we also write Sj as 5. 



4.2. Existence of a Weak Entropy Solution. For < s < t < Tq, define Q s j ■= 0,n{x E [s,t]}, 
S s = {l(~){x = s}, and = Df]{s < x < t}. By the definition of weak entropy solutions (Definition 
II. ip . a pair of bounded measurable functions (g, U) = (g(x), U(x,y)) is a weak entropy solution of 
problem (jl.5p provided that 
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• For any V 6 Qf 3 (M 2 ), 

Fg(U) := / (pud x ip + pvdyip) dy dx + / puipdy — / pmf)dy = 0; (4.5) 

• For any tp G C£°(R 2 ), 
:= y ((/ou 2 +p)d x ij) + puvdyip) dy dx + y" (pu 2 +p)ipdy 

(pu 2 + p)if; dy — p / V n i ds = 0; (4-6) 



• For any V G Cg°(K 2 ), 



/*([/) := / [puvdxij) + {pv 2 + p)d y ip) dy dx + \ puvipdy — j puvipdy 

-p j ifm 2 ds = 0; (4.7) 



I (pu(E + ^)8 X ^ + pv(E + -)d y if) dy dx 
Jn.,t P P 

+ pu{E+ -)ipdy- I pu{E + -)ipdy = 0; 
Jz s P JSt P 



• For any ip G C^° m2 ^ 



o 



Jt(U) := 



(4.8) 



• For any -0 G Cq^M 2 ) that is nonnegative, 

£*([/) := / (puSd x tp + pvSdyi)) dy dx + / pu^ dy - / dy < 0. (4.9) 

4.2.1. Estimate on the Total Strength of the Removed Fronts. For any approximate solution (g s , L^ 5 ), 
we set 

ft 5 := {x > 0,y > g 5 } 

and 

:= {y = /(x)}. 

For < s < t < T , define 

nf jt := ^n{xe [ s , t]}, sf = n = a}, r* it = r 5 n {x g [ Sj t}}. 

We note by our construction of approximate solutions that U s may not be a weak entropy solution 
of the Euler equations (jl.ip in f2 5 , since there are possible errors introduced by the approximating 
rarefaction wave via several fronts, and removing weak fronts of higher generations. In the following 
we will estimate these errors and show that they actually vanish as 5 — > 0. The analysis is again 
quite similar to [H]. We first list below Lemma 6.5 in [11] for later reference. 

Lemma 4.3. Let Q m denote the set of all fronts of generation m, and let T m denote the sum of 
the strengths of fronts of generation m: T m = J2 aj eg m \ a j\- Then T = ^2^ t =i%m an ^ 

T m < C{AKT) m 

for some constant C . In particular, for m = N + 1, we have T/v+i < C6 (cf. (|3. 13j) ). 
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4.2.2. Exact Riemann Solutions. For a given approximate solution (g s , II s ), suppose as before that 
the collision/reflection "times" are x = t\ < t<i < .... For a fixed interval [t,-,t,-+i], set si = r,-. 
We solve the following initial-free boundary problem with i = 1 (cf. (|2.ip ): 

' d x W(U) + d y H(U) = 0, x>s 4 , y>sOr), 

U = U 5 , x = Si, y > g{si) := g 5 (si), 

P = P, X > Si, on t := {y = g(x)}, 

v = g'u, x > Si, on T. 

bince the "initial data" U s (si, •) is piecewise constant, the solution (g\,U\) is obtained by solving 
the Riemann problems. It can be solved up to x = S2 when two waves interaction or reflection off 
the free boundary occurs (if S2 > Tj+%, we set S2 = Tj+i)- Then we solve (52,^2) from problem 
(|4.10p with i = 2 (note that the initial data is U s (s2, ■)), up to some S3. Repeat this process, we 
obtain 

(gi,Ui) in [si,s i+ i) 

with U^-Jsj, Sj + i) = [tj, Tj + i). We can then define (g,U) piecewise in x E [tj,Tj + i) by (g,U) = 
(gi,Ui) for x £ [si,s i+ i). 

4.2.3. Error of Splitting Rarefaction Waves. Let U s be the approximate solution obtained from 
problem (|4.10p in [sj,Sj+i], with the approximating parameter 5. This means that the rarefaction 
waves in U are separated into many discontinuities; while there is no front to be removed since 
each front in U is of generation one. Also, by our rule of splitting rarefaction waves, the lowermost 
state of U s is the same as U. This implies that the corresponding free boundaries are the same, 
and both U, U s are defined in the same domain. The analysis below is similar to [11]. We present 
details here to show the ideas there still work for our free boundary problem. 

Suppose that there is a rarefaction wave in U with the lower state Ui and upper state U r . Then 
this rarefaction wave is replaced by a step function U s . There also holds 

\U 5 (x,y)-U(x,y)\ < 0(5) 

by our splitting process (it is zero for the points not in rarefaction wave fan). We also want to find 
the error in the L 1 -space. To this end, we note that there are at most ^jT^p steps, and the width 
of each step is at most (x — Si) AA, with AA the difference of characteristic speeds of two adjacent 
approximate fronts of each step — it is less than 0(5) (cf. (|2.1ip ). Using the mean value theorem 
(since we know uniform L°° bounds of U and U s ), and summing up for all rarefaction wave fans 
across x, we find 

\W(U 5 )(x,y) -W(U)(x,y)\dy < C [ \U 5 (x, y) - U(x, y)\ dy 

y>g(x) Jy>g(x) 

< 0(5) \Ur ~ U?\\ x ~ s i\ < 0(5)T.Y.(U 5 )\x - Sl \ = 0(5)\x - Si \. (4.11) 

k 

We note here that ^ fc \Uj? — Uf\ is actually controlled by the total variation of the initial data by 
using the property of the Riemann solution. Similar inequality also holds when W(U) is replaced 
by H(U). 
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4.2.4. Error of the Removing Weak Fronts. We then compare U s and U s in x S [sj, Si+i]. We note 
that both U 5 and U s satisfy the same initial data. The only difference between them is that some 
fronts in U & of generation N + 1 are ignored to obtain U s . Note that, by the removing fronts of 
generation N + 1, we always keep the lowermost state the same as before. This means that the 
free boundary of U s is the same as , hence still to be y = g(x) = g s (x). Consequently, U s is 
different from II s in x G (si,Sj+i) only in a number of wedges emanating from the discontinuities 
in U s (si, •), and in each wedge, the difference is bounded by the strength of the removing fronts a 
that are of generation N + 1. We also note the width of each wedge is controlled by 0(x — Sj). By 
Lemma 14.31 we then find 

f \W(U 5 )(x,y) - W(U 5 )(x,y)\ dy < C [ \U 5 {x,y) - U s (x,y)\ dy 
<0(\x- Si \) \a\ <0(6)\x-Si\. (4.12) 

Similar inequality is also true for H(U). 

4.2.5. Total Error of Approximate Solutions. Since U is obtained by the exact Riemann solvers for 
x £ [si,Si + \], there must hold (with Q s j and £ s , and T Sjf in the integrals replaced by fif t , Sf, and 
T s s t respectively, since we have shown that the free boundary of U is the same as II s ): 

f:; +1 (u) = o, g^ 1 (u) = o, fsi +1 (u) = o, j s s ; +1 (ii) = o, (u) < o. 

From (|4.1ip and (|4.12p . we also obtain that, for any x S [sj, Sj+i], 

\W(U 5 )(x, y) - W(U)(x, y)\dy< 0(5)\x - Si \, (4.13) 

/ \H(U 5 )(x,y)-H(U)(x,y)\dy<0(5)\x- Si \. (4.14) 

Therefore, as an example, we find (note that the boundary term involving the pressure p j r ipni ds 
canceled because the boundary is the same): 

|G*: +1 (£/ 5 )| = |G** +1 (£/ 5 )-G*: +1 (£/)| 

^ [ Sl+1 I I [W 2 (U S ) - W 2 (U))d x <f> + {H 2 (U 5 ) - H 2 {U))d y <l>\ dydx 

+ [ \{W 2 (U 5 )-W 2 (U))<p\dy+ f \(W 2 (U s )-W 2 (U))<p\dy 

< M f 1+1 f \W 2 {U 5 )-W 2 {U)\dydx + M f \H 2 (U 5 ) - H 2 {U)\ dydx 

+M f \W 2 (U 5 ) - W 2 (U)\dy + M [ \W 2 (U s )-W 2 (U)\dy 



< 0(S)(s i+1 - Si y + 0(6)(s i+1 - Sl 
where M := ||0|| w i lOO . Then we find 

oo 

\G T T ] + \U 5 )\ < 0(S) ((*+i " + (*+i " < 0(<5)((r i+ i - r 3 f + (r j+1 - r 3 )). 



i=l 
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Thus it is clear that 

itfTrS 



Gl(U d )\ < 0(5)(\t - s\ z + \t- s\) for any < s < t < oo, 



and 



We now need to prove 



limG*(C/ 5 ) = 0. 

5—s-O 



YimG t s (U s ) = G t s (U). 
<5->0 



4.2.6. Verification of Weak Entropy Solutions. Set 

4>\x,9) = <f>(x,9 + g s {x)), Wi(x,0) = W 2 (U s )(x, 9 + g s (x)), H S 2 (x,9) = H 2 (U s )(x, 9 + g s (x)), 
where W$ = W 2 (U 5 ) and H 6 2 = H 2 (U 5 ). Then we have 

GKU 5 ) = If [Wid^ + Hldy^dydx 

Js Jy>g 6 (x) 

+ [ {wU)\ x=s dy- I (W^)\ x=t dy-p I (jmtds 

Jy>g 5 (s) Jy>9 S {t) Jr s st 

= 11 (W$(dJ s -de4> 5 (g 5 y) + H s 2 d e ft)d9dx 

Js Jr+ 

+ [ {Wi4>% =s d9- I (W^% =t d9 

Jr+ JR+ 



-p I ' <f>(x, g S(x)) — {9 J[ x) y/l + ((g s )'(xW dx. 

Js v 1 + (iff ) 0*0) 

By Lemma [4. 1|, we know that g s —> g uniformly for x G [s, t]. Since U s is uniformly bounded and 
converges to U under the metric of C([s,i]; L 1 (R + )), then W s and H s are also uniformly bounded 
and converges to W and H respectively in the topology of C([s,i]; L 1 (]R + )). From these facts, one 
can easily use the Lebesgue dominant convergence theorem to show (with <j) = <fi(x, 9 + g{x))) that, 
as 5 — > 0, 

/ {wU%=sd9- I [wU% =t d9^ I {W 2 h\x=sd9- I [W 2 $)\ x=t d9, 
Jr+ Jr+ Jr+ Jr+ 

f I {Wldj s + H 5 2 d e 4> 5 ) d9dx -> f I {W 2 dJ + H 2 d e 4>) d9dx. 

Js JR+ Js Jr+ 

Since {(g^)'} is uniformly bounded, we may assume that (g s )' — h in the weak* topology of 
L°°(]R + ). Since g 5 — > g uniformly in [s,t], we find that (g s )' — > g' in the sense of distributions. 
Thus, we must have h = g'. Therefore, as 4>(x, g s (x)) — > 4>(x, g(x)) uniformly in [s, t], and (g s )' — 1 g' 
in the weak* we have 



f (f>(x, /)(/)'(*) dx -> f 4>{x, g)g'(x) dx. 

J S J S 
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We also find 

[Wide4> s {g 5 )' - W 2 d e h') dOdx 



/'/ 

Js Jr- 



f [ Wl{de¥ - d e $) (g 6 )' dddx + [ [ Wid e 4>{(g 5 )' - g') dOdx 

Js JR+ Js Jr+ 

(W$ - W 2 )d e h'{x) dOdx. 



Using the boundedness of an d {(fl l<5 )'}> an d the uniform convergence dg(f> & — > d$<j>, the first 

integral in the right-hand side goes to zero as S — > 0. The third one converges to zero follows 
directly from W 2 — > W 2 in C([s, t]; L 1 (M. + )). For the second integral, it can be written as 



(Wi - W 2 )d e <P((g d )' - cf) dOdx + / / W 2 de<P{(g d Y - g') dOdx. 

Js Jr+ 

By the boundedness of (g 6 )' — g' , the first one then converges to zero; for the second one, using 
again (g s )' g' in the weak* topology of L°°. Then we have also proved 



W$do4> s (g 6 )' d6dx -»• f [ W 2 d e H' d9dx. 

Js JR+ 



Hence, we have 



limGKU 5 ) = I I {W 2 (dJ-d 4>g') + H 2 d $) dOdx 
<^o J s J R + 

+ [ {W 2 4>)\ x=s d0- [ (W 2 $)\ x=t d6-p f <t>(x,g{x))g'{x)dx 
Jr+ Jr+ Js 

= Gl(U) 

by a change of variables (x, y) = (x, 9 + g(x)). 

Therefore, we have proved G*({7) = as desired. Similarly, we can conclude 

Ej(U)<0, F t s (U) = 0, I t s (U) = 0, Ji(U) = 0, 

and hence the limit (g,U) obtained from the approximate solutions (g s ,U s ) is actually a weak 
entropy solution to problem (jl,5|) . 

It is clear that g should satisfy the estimate listed in Theorem ll.il as guaranteed by Lemma [4. 11 
To show || (Z7 — U_ + )(x, • ) 1 1 b v — ( ~' £ ' we no * e we nave P rove d 

\\(U 5 -U + )(x,-)\\ B v<Ce. 
Then, by Helly's theorem, without loss of generality, we may assume 

{U 6 - U + )(x, •) -»• (U - U + ){x, •) pointwise 
for some (U — U_ + )(x, •) so that 

\\(U -U + )(x,-)\\ By < Ce as 5^0. 
However, by uniqueness of the pointwise limit, we must have U = U. This completes the proof. 
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5. Asymptotic Behavior of Weak Entropy Solutions 

Finally we discuss the asymptotic behavior of the weak entropy solution (g, U) as x — > oo. 

For any given 5 > and the corresponding approximate solution (g s ,U s ), we know that there 
are a finite number of fronts and collisions/reflections. Thus, there exists x$ > so that, for x > x$, 
there is no collisions and reflections. Suppose then that there are m + 1 different states {L^}^_ 
from the upper to lower. It is obvious that Uq = U_q (cf. Remark II. ip . and there is mo with 
1 < mo <mso that each pair (Uj_ l , Uj) (j = 1, . . . , mo) is connected by a discontinuity of the 
first characteristic family, while each (Uj_ l , Ujj), j = mo + 1, . . . , m, is connected by a characteristic 
discontinuity (of the second and/or third characteristic family). Since no fronts interact, it is only 
possible that, for mo > 1, all the discontinuities of the first family must be rarefaction waves; for 
mo = 1, this discontinuity might be a shock or a rarefaction wave. For states Uj (j = mo+1, . . . , m), 
the pressure must be p by the boundary condition and the Rankine-Hugoniot jump conditions of 
characteristic discontinuities. 

Now we solve the free boundary Riemann problem of the Euler equations (jl.ip with the initial 
data U = U_q and the boundary condition p = p on the free boundary y = kx. Suppose that the 
solution is given by Uoo = ^4(a;U_ ). Then k = and the resulting 4-wave is a shock if p > p Q , 

and a rarefaction wave if p < p . It is clear that both (j = mo + 1, . . . , m) and (g ) should be 
7^ for all x > xs. 

We now define p s (£, 9) = + xg + l,0 + g s (£ + xs)) for £ > and 9 > 0. It is easy to see 
that \ps(£, -)l an d T.V.(p$)(£, •) = \p — p\ are bounded for all 5 > and given £. Thus, by Helly's 
theorem, there is a subsequence (still denoted as 5) so that lims->o ps(£,, 9) = p{9) = p for 9 > 
pointwise. This should imply that, for a.e. 9 > and the weak entropy solution U = (u,v,p,p), 

lim p(x, 9 + g(x)) = p. 

x — >oo — 

Similarly, we have 

lim — (x, 9 + g(x)) = 

x—i-oo u Uoo 

It then follows from q' = — that 

lim g'(x) = — — . 

t — ^ rv^ 1 1 
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